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ABSTRACT 

We pick out the 35 brightest galaxies from Goto's E+A galaxies catalogue which 
^^ , are selected from the Sloan Digital Sky Survey Data Release 5. As E+As have ex- 

perienced starburst recently and quenched it abruptly, they have been considered as 
post-starburst galaxies. The spectra of E+As are prominently characterized by the 
strong Balmer absorption lines but little [On] or Ha emission lines. In this work we 
O ' study the stellar populations of the sample galaxies by fitting their spectra using 

ULySS, which is a robust full spectrum fitting method. We fit each of the sample with 
^ . 1-population (a single stellar population-a SSP) and 3-population (3 SSPs) models, 

' separately. By 1-population fits, we obtain SSP-equivalent ages and metallicities which 

correspond to the 'luminosity-weighted' averages. By 3-population fits, we divide com- 
ponents into three groups in age (old stellar population-OSP, intermediate- age stellar 
population-ISP, and yotmg stellar population- YSP), and then get the optimal age, 
metallicity and population fractions in both mass and light for OSP, ISP and YSP. 
Q^ ' During the fits, both Pegase.HR/Elodie3.1 and Vazdekis/Miles are used as two inde- 

pendent population models. The two models result in generally consistent conclusions 
as follows: for all the sample galaxies, YSPs (^ IGyr) make important contributions 
r^ ' to the light. However, the dominant contributors to mass are OSPs. We also recon- 

^^ , struct the smoothing star formation histories (SFHs) by giving star formation rate 

^^ ■ (SFR) versus evolutionary age. All the sample galaxies have low SFRs in the inter- 

mediate evolutionary stage. Eleven of the thirty-five E-|-As have high SFRs in the 
early evolutionary stage. However, another 11 have SFRs that are high during the late 
evolutionary stage. This might be due to the recently happened but abruptly trun- 
. ^ cated starburst in such galaxies. In addition, we fit the E-f A sample and 34 randomly 

^ . selected elliptical galaxies with 2-population (2 SSPs) model, which could divide the 

galaxy components into two groups in age (old component and young component). 
We obtain the equivalent age of old components for each of the E-l-A sample and 
elliptical galaxies. By comparison, the old components of E-|-As are statistically much 
younger than those of ellipticals. From the standpoint of the stellar population age, 
this probably provides an evidence for the proposed evolutionary link from E-|-As to 
early- types (E/SOs). 

Key words: galaxies: evolution - galaxies: formation - galaxies: stellar content 



1 INTRODUCTION galaxies or an ensemble of K-type stars (Mgsirs, Fe527o and 

Ca3934 3468 absorption lines but no emission l ines) and of A- 

There exists such a class oi galaxy as has optical spec- , ,' , , „ , , .-,■ \ JTT^ — 1 ; — 71 

° , , . ,. , type stars (strong Balmer absorption lines) (IDressler et al.l 

trum characterized by strong Balmer absorption lines, but ^ Therefore, the so-called E-f A or K+A galaxy is only 

no or insignificant emission m On or Ha. These galaxies i -n ^- c ^ ^^ j^ j^ ^ -n ^- r 

, ,,^ , , r , 1,1, ^ classification of spectrum but not a cfassincation of mor- 

are calfed E-f- As or K-fAs , for their spectra took like , , „, ... „ , ,-, , , ,. ,. 

pnology. ifie existing of strong Bafmer absorption lines m 
a superposition of passive continuum spectra of efiptical ,, ,.,.,., ,, ■ i . , 

tfie spectra indicates the recently experienced starburst or 

the presence of young stars in these galaxies. However, no 
* Email: wdu@nao.cas.cn detection of emission line in [Oil] rules out the ongoing star 
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formation in these galaxies. So E+As are considered as post- 
starburst galaxies. 

It is proposed that E+As are good candidates in the 
middle stage of the transformation from late-types (star- 
forming, gas-rich, disk-dominated galaxies) to early-types 
(quiescent, gas-poor, spheroid-dominated galaxies). They lo- 
cate at such a crucial evolutionary stage that it is essential 
to study them for a better understanding of the galaxy evo- 
lution and for exploration of the origin of the red sequence 
of galaxies. So far, studies on E+As and the evolutionary 
link with early-type galaxies have been presented in several 
wayS; 

lYang et al.1 1|2004 l2006l . 120081 ') have presented 21 E-|-As 
from high-resolution HST ACS and WFPC2 images. Based 
on surface photometry and 2-dimensional model fit, they 
studied the morphologies, color profiles, new star clusters 
formed during the starburst, the location on the fundamen- 
tal plane (FP) of elliptical galaxies, and the evolution of the 
FP zero point. By disturbed morphologies and locations in 
the field, they suggested that galax;y-galaxy mergers and in- 
teractions are one mechanism which triggers off the E-l-A 
phase. The discovery of LINER spectra in the cores of some 
E-|-As in the follow-up spectroscopy gave confidence in the 
propos ed evolutionary lin k (E-|-As -^ E/SOs). 

In lPracv et al.l ((20091), they have presented 10 relatively 
bright (bj ~ 18.4) and nearby (z < 0.2) E-l-A galaxies 
from 2dFGRS imaging and spectroscopic observations. They 
made combined photometric and spectroscopic study for the 
10 E-|-As. In photometry, the structures, morphologies and 
spatial distributions of colours have been investigated. Mor- 
phologically the sample galaxies have taken on a look of 
tidal features, which indicated a high rate of interactions. In 
addition, the radial surface brightness profiles of the sample 
galaxies have enhanced the convincingness to the proposed 
evolutionary link (E-|-As — >■ E/SOs). In spectroscopy, stellar 
populations have been constructed by fitting spectra with 
template. Young stellar populations they have obtained for 
all the sample galaxies showed significant rotation, and this 
phenomenon commonly existed among overall populations 
of e arly-types. 

iFalkenberg et al.l (|2009l ) compared the spectra of E-l-A 
galaxies with the template of E and SO spectra, and this 
comparison indicated that early-type galaxies are succes- 
sors of E-l-A galaxies. They used the galaxy evolution code 
to trace the evolution of the Lick indices and colors of the 
E-l-A galaxies. In t erms of internal structures of galaxies, 
lYamauchi fc Gotcil (|2005l ) have studied the radial and 2D 
colour properties of 22 E-|-As from Sloan D igital Sky Sur- 
vey Data Release 2 (|Abazaiian et al.ll2004l ) and suggested 
that E-l-A galaxies are post-starburst galaxies which have 
undergone a centralized s t arbur st arising from interactions 
and mergers. IChoi et al.l l|2003) compared the fundamen- 
tal planes (FP) of E-l-A galaxies with FP of UV-excess 
early-type galaxies. Their work suggested that E-l-A galax- 
ies might undergo the transition from 'blue cloud' to 'red 
sequence', and eventually migrate to the red sequence early- 
type galaxies . 

More studies will rely on E-f As, so a large number 
of E-f A galaxies are needed from various telecopes and 
surveys. So far, there are already several E-l-A catalogues 
which are separately s elected from Sloan Digita l Sky Sur- 
vey Data Release 5 (|Adelman-McCarthv et al.1 bOOTl ) by 



iGotol (l2007l), the 2dF Galaxy Redshift Survey (2dFGRS) 
bv lBlake et al. 112004). the Las Ca mpanas Red shift Survey 
fLCRS) bv IZabludoff et al.l l|l996i ) and so on. iBlake et al.) 
(2004) has selected sample of low- redshift (z ~ 0.1) E-|-A 
galaxies from the 221,000 galaxy spectra in 2dFGRS. They 
adopted two different selection techniques, based on three 
Balmer absorption lines (H5, H7, H/?) or solely the H(5 line. 
These tw o met h ods y ield 56 and 243 E-f-As, respectively. 
iZabludoff et al.1 (|l996l ) identify 21 low-redshift (z ~ 0.05 - 
0.13) E-|-As from the 11,113 galaxies in the LCRS by using 
criteria that the average of equivalent widths of H/3, 7, 5 is 
greater than 5.5 A and the equivalent width of [On] emis- 
sion line is less than 2.5A. Our sample is taken from Goto's 
E-l-A catalogue, and it will be described in next section. 

As stellar populations of galaxies hold a wealth of infor- 
mation about the formation history and evolution, we shall 
study stellar populations of the E-l-A sample. In addition, 
we expect to explore new pieces of evidence in the aspect of 
stellar populations for the proposed evolutionary scenario 
of E-f As. At present, stellar populations can be widely con- 
structed by the stellar population synthesis technique which 
is based on fitting photometric indices and colours or spec- 
troscopic features with stellar population models. In this 
work, we choose 35 brightest E-fA galaxies from Goto's cat- 
alogue (|Gotcill20C)7l ). and th en use ULySS , which is a stellar 
population synthesis code (jKoleva et al.l l2009al ). to study 
the stellar populations and star formation histories. Com- 
paring with the previous works on E-f As , our work is based 
on spectroscopy instead of photometry since spectra con- 
tain much wealthier information than images. Our sample 
will be larger in number, and taken from SDSS spectroscopic 
catalogue, which has acknowledged higher quality spectra. 

In this paper, we describe our E-f A sample in § 2. We 
model the stellar populations for the sample and make con- 
clusions on dominance of stellar populations in § 3. In § 
4, we verify the reliability and robustness of our solutions. 
We present star formation rates (SFRs) to reconstruct the 
smoothing star formation histories (SFH) for the sample in 
§ 5. Additionally, in § 6 we provide a piece of evidence for 
the proposed evolutionary scenario of E-f As. We summarize 
the work in 5 7. 



2 DESCRIPTION OF THE E+A SAMPLE 

iGotd l|2007l ) has established a catalogue of 564 E-f A galax- 
ies identified from all galaxy s pectra of the SDSS DR5 
(JAdelman-McCarthv et al.llJOOTl '). These E-f As are selected 
by the following criteria: 

(l)spectroscopic signal-to-noise (S/N) ratio > 10 per 
pixel (in the continuum of the r-band wavelength range) 

(2)EW(H5)>5.0A, EW(Ha)> -3.0A and EW([Oii])> 
— 2.5A (EW is short for 'equivalent width' and EWs of these 
three lines are measured by ICotd l|2003h . Here absorption 
line has a positive sign.) 

(3)galaxies with redshift between 0.35 and 0.37 should 
be excluded due to the sky feature at 5577A after (1) and 

(2). 

By these criteria, the redshifts, absolute magnitudes 
an d diameters of the whole E-f A galaxies in this catalogue 
of ICotd (I2OO7I) have distribution histograms as shown in 
Figure 1 (the dashed lines). 
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We picked the 35 brightest E+A galaxies as our sample 
from Goto's catalogue. We chose this sample to ensure that 
the objects are sufficiently bright for further spectroscopic 
study in our work. These 35 galaxies are nearby and the 
redshifts are all within 0.14 (z<0.14). They have Petrosian 
magnitudes between 14 and 16 in r-band, span a redshift 
range from 0.0381 to 0.1322 (Figure 1, the solid line in the 
left histogram), and have high spectroscopic S/N ratio in 
the continuum of the r-band wavelength range (S/N > 30). 
In Table 1, we list the major properties of the sample. 

In Figure 1, we pr esent t he di stribution histograms of 
the full 564 E+As in iGotol (120071 ) and the selected sub- 
sample in terms of redshift, absolute magnitude in mag, 
diameter in kpc, and diameter in arcsec. For easy compar- 
isons, we plot the histograms of the full 564 (dashed lines) 
galaxies and 35 sample galaxies (solid lines) in the same 
picture in distinctive line styles. Particularly, we have trans- 
formed distributions of the 35 sample galaxies to the equiv- 
alent distributions if the number of 35 is enlarged to 564 
. This transformation is just for a more direct comparison 
with the distribution of the full 564 E-|-As. Considering our 
selection criteria (brightest, high S/N and then nearby), dis- 
tributions of the full 564 and 35 galaxies in redshift, absolute 
magnitude and diameter are generally consistent. Such con- 
sistency convinces us that the 35 E-l-A sample are generally 
good representatives of the full 564 E-|-As. 



3 CONSTRUCTING STELLAR POPULATIONS 
FOR OUR SAMPLE 

The SFHs are always likely irregular, since they may be both 
triggered and truncated by environmental e ffects, which 
have some r andomness( Maver et al.l l200ll : iHaines et al.l 
l2007l : lBosch et al. 2008). For example, E-l-A galaxies started 
violent starbursts recently and stopped starbu rsts abruptly 
l|Couchet al.lll987l : iDressler et akl] 19831 . Il992l l. For explor- 
ing SFHs of galaxies, it is essential to study the stellar pop- 
ulations of galaxies which can be constructed by analysing 
the galaxies spectra. The average properties of all the stel- 
lar populations in a galaxy can be simply derived by 1- 
population (a single stellar population-a SSP) fit for the 
galaxy spectrum, and the obtained average properties are 
commonly described by SSP-equivalent age and metallicity. 

In this section, we fit each of the sample with 1- 
population model to derive the equivalent properties of the 
representative stellar population for each of the sample. The 
principle of our analysis method is to compare an observed 
spectrum with a linear combination of some non- linear com- 
ponents (SSPs). We adopt two independent stellar popula- 
tion models in this work. One is the Pegase-HR/Elodie3.1 
population model, and the other is the Vazdekis/Miles pop- 
ulation model. 

The Pegase-HR/Elodie3.1 (PE) population model uses 
the Elodie 3.1 stell ar library l|Prugniel fc SoubiranI l200ll : 
IPrugniel et al.l [20071), which is the newest version of the 
Elodie library, containing 1962 spectra of 1388 stars from 
the spectroscopic arch ive of Observatoire de Haute-Provence 
l|Moultaka et all 120041 '). The wavelength coverage is 3900 ~ 
68OOA, and the spectral resolution (FWHM) is 0.55 A. The 
principle of the data reduction, flux calibration and deter- 
mination of the atmospheric parameters are described in 



IPruexiiel fc SoubiranI ll200lh. Then the library is operated 
by the Pegase-HR code l|Le Borgne et al.l |200J), which al- 
lows a choice of different physical ingredients (initial mass 
function-IMF, evolutionary track, SFR). The isochrones are 
solar-scaled at different values of the total met allicity Z. The 
evolutionary tracks are taken from Padova 94 (|Bertelli et al.l 
Il994h. and the SSPs are computed with Salpeter IMF 
dSalpeteiJ (J1955'); 0.1 Mq ^ M ^ 125 Mq). Ultimately, the 
PE model consists of 476 SSPs, covering 0.01 ~ 20 Gyr for 
age and -2.30 ~ 0.70 dex for [Fe/H]. 

The Vazdekis/Mil es (VM) popula tion model ad vances 
the pr evious model in IVazdekisI (Il999|) and Vazdck is et al 
(200s). It uses the Miles library (|Sanchez-Blazauez et al 
l200eiV which has a 2.3 A spectral resolution (FWHM) and 
3525 ~ 7500 A wavelength range. This library is believed 
to be bette r flux ca librated than any other. Padova 2000 
isochrones (jGirardi c t al. 2000|) are used, which has hotter 
red gi ant branch, an d the SSPs are computed with Salpeter 
IMF l|SalpeteiJll955l '). The generated VM model consists of 
276 SSPs, covering 0.1 - 17.7 Gyr for age and -1.7 ~ 0.2 
dex for [Fe/H]. 

W e accomplish our analysis by ULySS IjKoleva et al.l 
|2009a|), which is an open-source software package enabling 
full-spectrum fit for the study of stellar populations of galax- 
ies and star clusters. In ULySS, an observed spectrum is fit- 
ted with a model expressed as a linear combination of non- 
linear components. A component is a non-linear function of 
age, [Fe/H], [Mg/Fe], and wavelength. The ULySS method 
minimizes x^ value when matching an observed spectrum 
with models. It is worth mentioning that ULySS adopts the 
appoach of injecting the line spread function (LSF) into the 
model to match the resolution between the model and the 
ob servation. More deta ils about ULySS have been presented 
in lKolevaetHI (|2009al ). 

3.1 Pre-treating: Resolution matching 

To compare SSP model with an observed spectrum, we have 
to firstly match resolutions of the model and the observa- 
tion since they are different here. The spectral resolution 
is characterized by the instrumental broadening or the line 
spread function (LSF). LSF for spectra is equivalent to the 
point spread function (PSF) for images. In practice, LSF is 
not necessarily a Gaussian and varies with wavelength. The 
details on LSF have been described in iKoleva et al.l (120081 . 
l2009al ). To match the resolutions by ULySS, we need to de- 
termine the relative LSF between the model and the obser- 
vation, and then inject this relative LSF into the model. Our 
sample is selected from SDSS, so we just need to choose any 
spectrum from SDSS observations. Here we use the spec- 
trum of M67, which is a composite spectrum of stars from 
SDSS observations in M67 clusters and already contained 
in ULySS package, as representative of SDSS observations. 
ULySS calculate the relative LSF by comparing this com- 
posite spectrum of M67 with that extracted from the PE 
model. Then we inject this relative LSF to the PE model 
to generate the resolution-matched PE model by the LSF 
convolution function in ULySS package. Similarly, for the 
VM model, we do the same steps as above to generate the 
resolution-matched VM model. In the following text, we will 
use the resolution-matched PE or VM model, although they 
are still named as PE model and VM model. 
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Table 1. List of the 35 E+A sample selected from Goto's catalogue. Column 1 is the Galaxy Number (GN) we named for 
each of the sample. Column 2 is the SDSS ID. The rest of columns represent Petrosian magnitude in r-band, S/N ratio in 
r-band, redshift, right ascension (RA) , decl i nation (DEC) and EW of H<5, [On] and Ha, respectively. It should be noted 



that EW values here are measured bv lCotd 1 12003) . 



GN 



lAU 
designation 



Petro 

mag(r) 



S/N 
(r) 



RA 



DEC 



EW 
H5 



EW 

[On] 



EW 
Ha 



1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 
35 



SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 
SDSS 



J130525, 
J162702, 
J123936, 
J120523, 
J115821, 
J145455, 
J155435, 
J080218, 
J141419, 
J084918, 
J120419, 
J105220, 
J161330, 
J210258, 
J100743, 
J233712, 
J125820, 
J153016, 
J124204, 
J110540, 
J113040, 
J104220, 
J080925, 
J133757, 
J115837, 
J133323, 
J083454, 
J104204, 
J105743, 
J104215, 
J162938, 
J092006, 
J134802, 
J124534, 
J135030, 



82+533530.3 
55+432833.9 
05+122619.9 
24+643029.6 
65+363820.4 
44+453126.5 
53+291319.9 
39+323207.8 
29-031111.5 
88+462252.5 
07-001855.8 
44+054941.5 
18+510335.5 
87+103300.6 
62+554934.5 
77-105800.4 
71+613039.4 
06+373346.1 
66+150905.6 
70+055954.2 
84+562910.0 
14+564855.7 
07+305652.8 
98+654410.5 
72-021710.9 
78+532130.1 
79+251258.5 
13+293323.4 
93+123539.9 
17+140636.3 
26+302924.9 
43+015807.7 
18+020405.6 
16+402559.1 
77+012804.7 



14.49 
14.64 
14.69 
14.93 
14.97 
15.04 
15.22 
15.22 
15.24 
15.30 
15.30 
15.34 
15.38 
15.42 
15.43 
15.43 
15.49 
15.57 
15.59 
15.64 
15.65 
15.68 
15.70 
15.76 
15.79 
15.80 
15.82 
15.85 
15.91 
15.96 
15.96 
15.97 
15.97 
15.98 
15.99 



50.4 
45.0 
43.4 
41.4 
48.3 
40.6 
45.7 
45.1 
38.9 
47.4 
40.1 
46.1 
33.9 
49.9 
54.5 
44.4 
40.6 
53.1 
43.9 
49.8 
52.9 
39.5 
44.7 
42.3 
37.6 
33.4 
38.8 
42.7 
43.0 
53.9 
36.8 
50.5 
38.6 
43.0 
39.1 



0.0381 
0.0463 
0.0409 
0.0822 
0.0656 
0.0366 
0.0943 
0.0382 
0.0471 
0.0410 
0.0938 
0.0412 
0.0336 
0.0928 
0.0452 
0.0783 
0.0905 
0.0775 
0.0872 
0.0543 
0.0622 
0.0909 
0.0836 
0.0665 
0.0880 
0.1023 
0.1322 
0.0398 
0.1195 
0.0543 
0.0971 
0.0850 
0.0678 
0.0822 
0.0728 



196.35761 
246.76067 
189.90021 
181.34687 
179.59021 
223.73102 
238.64808 
120.57664 
213.58040 
132.32870 
181.07948 
163.08521 
243.37578 
315.74529 
151.93174 
354.30318 
194.58634 
232.56693 
190.51944 
166.41961 
172.67019 
160.58393 
122.35451 
204.49158 
179.65718 
203.34913 
128.72832 
160.51723 
164.43308 
160.56321 
247.40944 
140.02682 
207.00909 
191.39237 
207.62819 



53.591759 
43.476095 
12.438888 
64.508254 
36.639018 
45.524055 
29.222199 
32.535504 

-3.1865321 
46.381272 

-0.31550661 
5.8282184 
51.059881 
10.550177 
55.826261 

-10.966755 
61.510986 
37.562823 
15.151574 
5.9984044 
56.486122 
56.815500 
30.948015 
65.736242 

-2.2863741 
53.358373 
25.216253 
29.556503 
12.594419 
14.110087 
30.490256 
1.9687923 
2.0682588 
40.433093 
1.4679767 



5.06 
5.68 
5.76 
6.79 
5.71 
5.16 
5.89 
5.54 
5.82 
8.34 
7.98 
6.21 
7.57 
5.19 
5.91 
6.66 
7.28 
7.10 
7.40 
6.89 
5.60 
5.50 
5.48 
8.02 
6.08 
6.16 
6.09 
6.44 
8.36 
6.33 
5.72 
6.01 
6.66 
6.18 
7.12 



1.12 
-1.42 
0.86 
0.90 
1.41 
0.96 
-0.60 
-1.75 
-1.12 
-1.16 
-1.18 
-0.63 
1.06 
0.30 
-2.06 
0.05 
-1.15 
-1.27 
0.63 
-1.42 
-0.98 
-1.78 
-1.86 
-1.36 
0.68 
-0.42 
1.04 
-1.03 
1.29 
-0.46 
-0.69 
-0.54 
-0.02 
0.80 
-1.70 



2.38 
0.68 
1.46 
1.78 
2.49 
1.00 
1.29 
1.54 
0.32 
1.91 
2.09 
1.92 
0.55 
1.81 
1.58 
1.76 
-1.81 
-2.82 
2.10 
1.01 
2.16 
1.18 
0.83 
2.29 
1.75 
1.06 
1.47 
1.52 
2.48 
2.16 
1.51 
-1.52 
2.11 
1.70 
0.63 



200 






150 




■ 


Number 




1 


■ 


50 





H 


-T -L, 



250 
200 

150 

100 

50 




F 


— 


^Z" 


n 



-0.1 0.0 0.1 0.2 0.3 
redshift 



-24 



-22 -20 -18 



200: 

150; 

100; 
50; 

0: 



absolute magnitude, mag 



r- 


- 


■ — 


> 


- -Li 



10 20 30 40 
diameter, kpc 



400 
300 
200 

100: 

o; 



■J L. 

1 rr 


-n .n 


: 



10 15 20 



diameter, arcsec 



Figure 1. Comparison of distribution histograms of some major properties of E+As. We compared the redshift, absolute magnitude in 
mag, diameters in kpc, and in arcsec for the 35 E+A sample (solid lines) and the full 564 E+A galaxies (dashed lines). We use bin=0.02, 
1.0 mag, 5.0 kpc, and 2.0 arcsec, respectively. For the 35 E+A sample, the mean values of redshifts, absolute magnitudes in mag, and 
diameters in kpc and in arcsec are respectively 0.0705, -21.86 mag, 14.27 kpc and 11.27 arcsec. For the full 564 E+A galaxies, they are 
respectively 0.1414, -21.73 mag, 15.18 kpc and 3.34 arcsec. 
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Table 2. Comparison of the mean SSP-equivalent age (Mage), 
[Fe/H] {MIp^/JJ^) and velocity dispersion (Ma) for the 35 bright- 
est E+As and the 564 E+As, illustrated respectively from the PE 
and VM model . 



Sample 



35 E+As 
564 E+As 



Mage(Gyr) 
PE VM 



0.87 
1.01 



1.14 
1.28 



^^Fe/HlCdex) 
PE VM 



-0.04 
-0.11 



-0.38 
-0.38 



Ma(km/s) 
PE VM 



159 
144 



141 
129 



3.2 A single stellar population for our sample 

The first step towards reconstructing the SFH of an ob- 
ject is to calculate the SSP-equivalent properties by fit- 
ting the spectrum with a SSP. This derived SSP represents 
the 'luminosity- weighted' epoch of star formation. There- 
fore, the SSP-equivalent properties are corresponding to 
the 'luminosity-weighted' average over the distributions and 
provide a general view of the stellar populations of the 
galaxy. 

In this section, we make a SSP model fit for each galaxy 
by using ULySS. Both the PE and VM model are indepen- 
dently used as population models. Finally, we obtain the 
SSP-equivalent age and metallicity for each galaxy. 

We show the results of a SSP fit with the PE model 
in Table 3 and with the VM model in Table 4. In both 
tables, the errors on the parameters are the formal 1- 
a errors, computed from the covariance matrix by the 
ULySS calling function: MPFIT which is a widely used 
fitting function. MPFIT gives out the optimal parameters 
from the best fit and the 1-cr errors on them. More de- 
tails on the errors can be obtained from the MPFIT algo- 
rithm (ihttp://cow. physics.wisc.edu/~craigm/idl/idl. html ). 
For a comparison, we show distribution histograms of SSP- 
equivalent age, metallicity and velocity dispersion of the 
sample in Figure 2 (solid lines) . Figure 2 shows that no mat- 
ter which model is used for fit, the derived SSP-equivalent 
ages are concentrating on about 1 Gyr, the metallicities dis- 
tribute mostly between -0.50 and +0.20 dex, and the velocity 
dispersions cover from 100 to 220 km/s. Such results might 
indicate that the populations aged around IGyr would be 
the 'light-weighted' components in E+A galaxies. 

Then we make a SSP fit for the total 564 E+A galaxies 
from Goto's catalogue with both the PE and VM model, 
and we also present the distribution histograms of the re- 
sults in Figure 2 (dot-dashed lines). We calculate the mean 
SSP-equivalent age, metallicity and velocity dispersion and 
tabulate them in Table 2 for clear comparisons. 

The results from both models are not exactly the same. 
The differences might arise from the difference and indepen- 
dence of the PE and VM model. The SSPs of the two models 
are generated by different population synthesis codes: Pe- 
gase.HR and Vazdekis, different stellar libraries: EL0DIE3.1 
and Miles, and different choices of the physical ingredients: 
evolutionary tracks, etc. However, the comparison in Figure 
2 and Table 2 shows that the results from the PE and VM 
model are generally consistent. 

In addition, both Figure 2 and Table 2 show the con- 
sistency between our 35 E+A sample and the full 564 E+A 
galax;ies in the aspect of the SSP-equivalent age, metallicity 
and velocity dispersion. This consistency further proves that 
our 35 E+A sample is a good representative of the full 564 
E+A galaxies. 



Table 3. Results of a SSP fit with the PE model. Column 1 is 
the Galaxy Number (GN). Column 2 is the x^ value from the 
best fit. Column 3 and 4 are respectively the SSP-equivalent age 
in Gyr, [Fe/H] in dex and the errors on them. 



GN X SSP-equivalent age 

(Gyr) 



SSP-equivalent [Fe/H] 
(dex) 



1 


0.43 


1.10+0.03 


0.02+0.03 


2 


0.64 


1.12+0.04 


-0.04+0.05 


3 


0.58 


1.09+0.04 


-0.02+0.05 


4 


0.67 


0.82+0.03 


0.16+0.04 


5 


0.44 


0.80+0.02 


0.12+0.03 


6 


0.61 


0.86+0.02 


0.17+0.03 


7 


0.59 


1.03+0.03 


-0.05+0.05 


8 


0.51 


0.93+0.02 


0.14+0.03 


9 


0.65 


1.15+0.06 


-0.10+0.05 


10 


0.53 


0.67+0.02 


-0.24+0.04 


11 


0.69 


0.71+0.02 


-0.19+0.06 


12 


0.54 


0.99+0.03 


-0.07+0.04 


13 


0.66 


0.06+0.00 


-0.07+0.06 


14 


0.61 


1.17+0.05 


0.02+0.03 


15 


0.67 


1.00+0.03 


-0.12+0.04 


16 


0.53 


0.80+0.02 


0.09+0.04 


17 


0.69 


0.82+0.03 


-0.68+0.08 


18 


0.56 


0.74+0.02 


-0.06+0.04 


19 


0.55 


0.73+0.02 


0.03+0.04 


20 


0.53 


0.78+0.02 


-0.08+0.05 


21 


0.64 


0.87+0.02 


0.08+0.03 


22 


0.64 


1.13+0.05 


-0.07+0.05 


23 


0.58 


1.00+0.03 


-0.01+0.05 


24 


0.62 


0.71+0.02 


-0.11+0.05 


25 


0.65 


0.88+0.03 


0.06+0.04 


26 


0.72 


0.88+0.04 


-0.01+0.06 


27 


0.64 


0.91+0.03 


-0.02+0.06 


28 


0.59 


0.94+0.03 


-0.05+0.04 


29 


0.57 


0.67+0.02 


-0.20+0.05 


30 


0.65 


0.87+0.02 


-0.14+0.03 


31 


0.74 


1.06+0.05 


-0.06+0.06 


32 


0.91 


0.96+0.04 


-0.14+0.05 


33 


0.64 


0.83+0.03 


-0.09+0.05 


34 


0.62 


0.82+0.02 


0.11+0.04 


35 


0.63 


0.68+0.02 


0.12+0.04 



3.3 The Composite Stellar Population for our 
sample 

It is too simple to fit the spectra of galaxies only using a 
SSP, since the galaxies have generally a complex SFH. Re- 
tracing the SFR along the history is a fundamental piece 
of information to understand the physics of the galaxies. In 
principle, one can obtain such information by directly fit- 
ting the galaxy with a positive linear combination of many 
SSPs, which means a composite stellar population (CSP) fit 
for the galaxy, but such an approach would be unstable be- 
cause of the degeneracies between the SSP components. To 
circumvent these degeneracies, ULySS takes effective mea- 
sures. It starts with the simple physical assumptions, such 
as the presence of an old stellar population, then divide the 
time axis into intervals by setting limits in two or more in- 
tervals. More details o n dealing with degene racies by ULySS 
have been analysed in Kolev a et al.l l|2009ah . 

In our work, we analyse the spectrum of each galaxy 
in terms of three epochs which means to fit each spectrum 
with 3 stellar populations: an old stellar population (OSP), 
an intermediate-age stellar population (ISP), and a young 
stellar population (YSP). They are respectively chosen from 
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Figure 2. Parameters distributions from a SSP fit with the PE model and the VM model. The solid lines show the distribution histograms 
of the 35 E+A sample, and the dot-dashed lines are for the full 564 E-|-A galaxies. We present the SSP-equivalent age, metallicity and 
velocity dispersion distribution of them. The top panel shows the results from the PE model, and the bottom is from the VM model. 
In the top panel, for the 35 E+A sample, the mean SSP-equivalent age, [Fe/H] and velocity dispersion are 0.87 Gyr, -0.04 dex and 159 
km/s; for the 564 E-l-A galaxies, the mean values are respectively 1.01 Gyr, -0.11 dex and 144 km/s. In the bottom panel, for the 35 E-l-A 
sample, the mean values are respectively 1.14 Gyr, -0.38 dex and 141 km/s; for the 564 E-l-A galaxies, the mean values are respectively 
1.28 Gyr, -0.38 dex and 129 km/s. These mean values are tabulated in Table 2 for a direct comparison. 



the following 3 boxes divided by different limits of age and 
metallicity. 

When we are using the PE model, we choose the boxes 
as follows: 

(i) For the OSP box, we fix its age to 12Gyr and set the 
[Fe/H] free (we fix the age to maintain the degree of freedom 
at the minimum, and meanwhile to see whether such an old 
stellar component can be discovered in E-fAs). 

(ii) For the ISP box, we make the age span a range 
from 1 to 5 Gyr and the [Fe/H] span from -1.00 to 0.20 dex 
(for several galaxies, the ISP component is pushing on such 
[Fe/H] limits, so we set the [Fe/H] free for these galaxies). 

(iii) For the YSP box, wo make the age span a range 
from 0.01 to 1 Gyr and the [Fe/H] span from -1.00 to 0.70 
dex. 

We fit each of the sample with 3 populations (OSP + 
ISP -I- YSP) by using ULySS, which gives the optimal results 
from the best fit by adjusting the set of parameters. It gives 
the optimal age, [Fe/H], light fraction (LF), mass fraction 
(MF) and errors on them for each component. During the 
fit, the boxes are suitable for all the sample galaxies except 
for several ones which have their optimal paramters reaching 
to the boundaries of boxes. For such cases, we only give their 
upper or lower limits. 

Results from this 3-population (3 SSPs) fit are gener- 
ally better than those from 1-population fit (§ 3.2) and 2- 
population (2 SSPs) fit (In § 6, we will additionally make 
the 2-SSP fit for each of the sample to compare EH- As with 
ellipticals) . The x^s are lower and the residuals, especially at 
the absorption features, are smaller. We show the exemplary 
fit and the residuals of Galaxy number 1 as a representative 
in Figure 3. Table 5 presents the optimal parameters and er- 



rors on them from this 3-SSP fit with the PE model, and the 
optimal age, [Fe/H], LF and MF are tabulated for each com- 
ponent. The errors on the optimal age and [Fe/H] in Table 
5 are obtained from the MPFIT fitting process similarly as 
mentioned for errors in Table 3, and errors on the light frac- 
tion and mass fraction are determined by a multiple linear 
regression routine called by ULySS. 

In Table 5, we find that 12Gyr-old components can be 
detected in all of the 35 galaxies. Very young(<0.13 Gyr) 
components are discoverd in 3 of the sample (Galaxy num- 
ber 3, 13 and 17), which indicates a quite recent starburst 
in the three E-l-A galaxies. For the majority of the sample 
(23 out of 35), YSPs dominate the light, since they have 
much larger light fraction than ISPs and OSPs. For 10 of 
the remaining 12 galaxies, the ISPs dominate the light. For 
the mass, all the 35 galaxies except 2 are dominated by the 
OSPs, since they have much higher mass fraction. To suma- 
rize, most of E-|-As are dominated by YSPs in light and some 
are dominated by ISPs in light. While the OSPs are not the 
dominant contributors to fight, but they are the dominant 
contributors to mass, statistically. This conclusion can be 
obtained more clearly from Figure 4 and Figure 5, which 
will be described later. 

For comparison we repeat fitting each of the sample 
with 3 SSPs from the VM model. To adjust to the narrower 
coverage of the VM model, we change the limits for boxes 
as follows: 

(i) For the OSP box, we still fix its age to 12Gyr and 
set the [Fe/H] free (we fix the age to maintain the degree 
of freedom at the minimum, and meanwhile to see whether 
such an old stellar component can be discovered in E-fAs). 

(ii) For the ISP box, we make the age span a range from 
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Table 4. Results of a SSP fit with the VM model. Column 1 is 
the Galaxy Number (GN). Column 2 is the x value from the 
best fit. Column 3 and 4 are respectively the SSP-equivalent age 
in Gyr, [Fe/H] in dex and the errors on them. 



GN 


2 
X 


SSP-equivalent age 


SSP-equivalent [Fe/H] 






(Gyr) 


(dex) 


1 


0.50 


1.39±0.05 


-0.23±0.01 


2 


0.63 


1.38±0.05 


-0.37±0.02 


3 


0.64 


1.39±0.06 


-0.38±0.02 


4 


0.76 


1.14±0.02 


-0.33±0.03 


5 


0.53 


1.02±0.02 


-0.26±0.03 


6 


0.67 


1.38±0.05 


-0.35±0.02 


7 


0.68 


1.38±0.05 


-0.35±0.01 


8 


0.59 


1.14±0.01 


-0.15±0.02 


9 


0.69 


1.38±0.07 


-0.29±0.02 


10 


0.57 


0.75±0.01 


-0.34±0.03 


11 


0.70 


0.84±0.02 


-0.38±0.03 


12 


0.57 


1.14±0.01 


-0.28±0.02 


13 


0.77 


0.60±0.01 


-1.33±0.06 


14 


0.70 


1.39±0.05 


-0.26±0.01 


15 


0.72 


0.96±0.01 


-0.05±0.02 


16 


0.60 


0.90±0.02 


-0.10±0.02 


17 


0.72 


1.49±0.02 


-1.21±0.03 


18 


0.73 


1.19±0.01 


-0.83±0.03 


19 


0.64 


0.89±0.01 


-0.28±0.02 


20 


0.62 


0.84±0.02 


-0.11±0.03 


21 


0.73 


0.86±0.02 


0.12±0.02 


22 


0.69 


1.38±0.06 


-0.33±0.02 


23 


0.62 


1.38±0.04 


-0.38±0.01 


24 


0.63 


O.SOitO.Ol 


-0.21±0.02 


25 


0.76 


1.17±0.02 


-0.38±0.03 


26 


0.87 


1.15±0.02 


-0.50±0.03 


27 


0.74 


1.60±0.18 


-0.67±0.02 


28 


0.63 


1.16±0.01 


-0.35±0.03 


29 


0.66 


0.75±0.01 


-0.27±0.02 


30 


0.69 


0.96±0.02 


-0.31±0.02 


31 


0.70 


1.37±0.06 


-0.39±0.02 


32 


0.96 


1.87±0.06 


-0.73±0.03 


33 


0.73 


0.94±0.02 


-0.23±0.03 


34 


0.66 


1.14±0.01 


-0.38±0.03 


35 


0.67 


0.88±0.02 


-0.30±0.03 



1 to 5 Gyr and the [Fe/H] span from -1.00 to 0.10 dex (Here 
we choose 0.10 dex because of the small upper limit (0.20 
dex) of the VM model for [Fe/H]; for several galaxies, we set 
the [Fe/H] free since the ISP component is pushing on the 
[Fe/H] limits). 

(iii) For the YSP box, we make the age span a range 
from 0.1 to 1 Gyr and the [Fe/H] span from -1.00 to 0.20 dex 
(no younger populations than O.lGyr and no richer popula- 
tions than 0.20 dex are contained in the VM model). 

In Table 6, we list the results for each component from 
the 3-population fit with the VM model. All of the sample 
galaxies have 12Gyr-old populations. Similarly as in Table 
5, very young populations (< 0.13) are detected in the same 
three galaxies (Galaxy number 3, 13 and 17). All of the 
results and conclusions from 3-population fit with the VM 
model (Table 6) have a good agreement with those from the 
PE model (Table 5). This general agreement could also be 
seen from Figure 4 and Figure 5, which will be described 
later. 

For a clear comparison of the results of the 3-population 
fit with the two independent population models, we show 
the distribution histograms of the light fraction (Figure 4) 



and mass fraction (Figure 5) of the YSP, ISP and OSP for 
all of the 35 E-fA sample. In both figures, the top panel 
shows results from the PE model, and the bottom panel 
shows results from the VM model. Light fraction distribu- 
tions of the YSP, ISP and OSP are presented from left to 
right, respectively. Obviously, we can conclude that no mat- 
ter which model is used for fit, in a statistical sense, YSPs 
contribute the most to the light, then the ISPs, and finally 
the OSPs; however, OSPs contribute the most to the mass, 
then the ISPs, and finally the YSPs. This conclusion can 
be obtained more clearly from Table 7, which presents the 
mean light contributions, mass contributions and ages of 
YSP, ISP and OSP by using the PE and VM models. YSPs 
contribute about 50% of the light, and ISPs and OSPs con- 
tribute another 50%. However, OSPs contribute about 70% 
of the mass, and ISPs and OSPs are presented in the rest 
30%. For the mean age, YSP of the sample is concentrat- 
ing on a mean age of 0.54 Gyr with the PE model and 0.46 
Gyr with the VM model, the ISPs are concentrating on the 
mean age of 1.82 Gyr with the PE model and 1.82 Gyr with 
the VM model. The two independent models give generally 
consistent results. 



4 RELIABILITY OF THE SOLUTIONS 

It is legitimate to question the reliability of the solutions be- 
cause of degeneracies, so in this section we are demonstrating 
the reliability of our solutions in terms of the robustness of 
ULySS techqiue, the validity of our solutions, the sensitivity 
to the population model and the comparison with another 
work. 



4.1 Robustness of synthesis technique 

ULySS is a package to fit spectroscopic observations against 
a linear combination of non-linear model components con- 
volved with a parametric line-of-sight velocity dispersion 
(LSF function in § 3.1). It fits the full spectrum with popula- 
tion models to derive the physical properties of stellar pop- 
ulations. This method appears as an optimized alternative 
to methods based on spectrophotometric indices or spectral 
energy distribution (SED), as it ensures a high utilization 
ratio of all the useful information contained in the signal by 
fitting all the pixels. The algorithm has been optimized in 
the mathematical details to improve the precision, robust- 
ness and performance, which finds out the optimal values of 
unknown parameters and the errors on them. The optimal 
age and metallicity for each component result from the MP- 
FIT function performed by the Levenberg-Marquardt least- 
squares minimization, and the optimal weight and luminos- 
ity weight of each component and the errors on them are 
determined by a multiple linear regression fit between the 
components and the observations (errors on the optimal age, 
[Fe/H], mass fration and light fraction of the components are 
all given in Table 5 and 6). It is convinced that results from 
ULySS method are consistent wit h those from only fitting 
Lick indices, but are more precise (jKoleva et al.l l2008l. 

It is an old problem that there are multiple solutions 
existing in stellar population synthesis. The multiplicity is 
caused by a combination of three factors: algebraic degen- 
eracy, intrinsic degeneracies of stellar populations, and the 
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Figure 3. Best fit witli tiie PE model and the residual spectrum for Galaxy number 1. The top panel shows the spectrum in black and 
the best fit in blue (both are almost superimposed and the black line can be seen only when zooming on the figure), the cyan line is 
the multiplicative polynomial to absorb the effects of an un precise flux calibrati on and of the Galactic extinction. The importance and 
details about the multiplicative polynomial are discussed in lKoleva et al.l | |200S|) . The red regions are rejected from the fit (rejection of 
flagged telluric lines and automatic rejection of outliers). The bottom panel is the residuals from the best fit. The continuous green lines 
mark the l-cr deviation, and the dashed green line is the zero-axis. 
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Figure 4. Distributions of the 3 populations with respect to light fraction (LF). The top panel shows LF distributions of YSP (left), 
ISP (middle), and OSP (right) from 3-population fit with the PE model. The bottom panel shows the same but with the VM model. For 
the top panel, the mean LFs of all YSPs, ISPs and OSPs are 51%, 30% and 19%. For the bottom panel, the mean LFs are respectively 
46%, 35% and 19%. These values are tabulated in Table 7 for a direct comparison. 
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Figure 5. Distributions of the 3 populations with respect to mass fraction (ME). The orientation is the same as that in Figure 9. For 
the top panel, the mean LF of all YSPs, ISPs and OSPs are 11%, 14% and 75%. For the bottom panel, the mean LFs are respectively 
12%, 25%, and 63%. These values are tabulated in Table 7. 
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Table 5. Parameters from 3-population (OSP + ISP + YSP) fit with the PE model for each of the sample. This table presents the 
population histories of the 35 E+As, reconstructed from 3 episodes of star formation. Column 1 is the Galaxy Number (GN) for the 
E+A sample listed in Table 1. Column 2 is x^ value from the best fit. The following are age (in Gyr), [Fc/H] (in dex), mass fraction 
(MF, in %), and light fraction (LF, in %) for each of the component (OSP, ISP, and YSP). The last line lists the results of the combined 
spectrum, which will be described later in § 4.2. 



GN x' 


OSP 


(fixed to 12Gyr) 


ISP 


YSP 




[Fe/H] 


MF 


LF 


Age [Fe/H] 


MF 


LF 


Age [Fe/H] MF LF 




(dex) 


(%) 


(%) 


(Gyr) (dex) 


(%) 


(%) 


(Gyr) (dex) (%) (%) 


1 0.40 


-0.22±0.64 


61.16±3.28 


16.46±0.88 


1.05±0.29-0.50±0.37 


14.69±0.24 


43.65±0.72 


0.99±0.33 0.52±0.66 24.15±0.95 38.89±1.57 


2 0.58 


0.21±0.11 


92.51±1.54 


38.80±0.64 


1.47±0.32-0.69±0.25 


4.72±0.32 


28.53±1.91 


0.51±0.19-0.38±0.28 2.77±0.11 32.67±1.28 


3 0.51 


0.25±0.08 


94.79±0.57 


42.27±0.25 


0.73±0.18-0.38±0.12 


4.95±0.08 


48.26±0.82 


0.09±0.09 0.30±0.77 0.26±0.02 9.47±0.59 


4 0.63 


0.26±0.18 


89.40±0.52 


28.11±0.16 


1.21±0.76-1.00±0.63 


2.54±0.39 


15.25±2.35 


0.58±0.10 0.06±0.23 8.06±0.32 56.64±2.24 


5 0.41 


0.21±0.61 


56.91±5.51 


9.38±0.91 


1.49±0.47 0.30±0.20 


28.04±1.32 


32.81±1.54 


0.60±0.07 -0.23±0.10 15.05±0.17 57.81±0.65 


6 0.57 


0.22±0.25 


61.76±6.75 


12.08±1.32 


1.48±0.27 0.41±0.14 


24.80±1.65 


32.67±2.17 


0.70±0.07-0.25±0.13 13.45±0.21 55.25±0.87 


7 0.57 


0.29±0.29 


80.51±0.97 


19.27±0.23 


1.55±0.27-0.84±0.36 


11.69±0.54 


42.19±1.95 


0.56±0.07 ^ 0.30 7.80±0.36 38.54±1.76 


8 0.50 


-0.06±0.19 


63.45±1.01 


13.19±0.21 


1.42±0.37-0.70±0.59 


8.89±0.78 


21.14±1.86 


0.73±0.04 ^ 0.30 27.66±0.83 65.67±1.96 


9 0.64 


0.36±0.22 


80.26±3.28 


19.23±0.79 


1.27±0.12-0.71±0.47 


14.47±0.49 


64.66±2.18 


^1.00 0.30±0.52 5.27±0.95 16.12±2.92 


10 0.47 


-0.19±0.45 


76.20±7.79 


21.30±2.18 


6.66±8.75 0.42±0.93 


13.34±6.60 


3.87±1.91 


0.33±0.03-0.14±0.08 10.46±0.04 74.84±0.29 


11 0.61 


0.13±0.43 


74.10±6.31 


18.66±1.59 


3.89±1.62-0.12±0.71 


17.80±2.68 


14.25±2.15 


0.34±0.05-0.10±0.15 8.11±0.07 67.08±0.58 


12 0.52 


0.25±0.45 49.78±12.80 


7.45±1.92 


1.84±0.94 0.52±0.14 


18.80±3.97 


13.53±2.86 


0.96±0.06 -0.31±0.10 31.42±0.38 79.02±0.96 


13 0.66 


-0.69±1.14 


51.97±9.16 


8.88±1.57 


7.87±5.97 0.69±1.17 37.90±12.95 


4.67±1.60 


0.08±0.01 -0.20±0.09 10.12±0.02 86.45±0.17 


14 0.55 


0.15±0.16 


79.83±3.15 


26.56±1.05 


1.38±0.29 0.49±0.12 


13.90±0.71 


30.04±1.54 


0.75±0.10-0.59±0.20 6.27±0.07 43.41±0.51 


15 0.64 


0.54±0.11 


69.86±2.05 


10.43±0.31 


1.28±0.24-0.55±0.46 


17.65±1.05 


50.10±2.98 


0.74±0.34 0.04±0.30 12.49±0.86 39.47±2.71 


16 0.50 


0.22±0.17 


87.77±1.30 


25.71±0.38 


^1.00 -0.48±1.53 


1.93±1.11 


9.61±5.52 


0.64±0.33-0.14±0.40 10.30±0.82 64.68±5.16 


17 0.57 


^-2.30 


18.09±1.40 


12.88±0.30 


^5.00 0.44±0.12 


78.98±0.32 


19.45±1.14 


0.12±0.01-0.95±0.22 4.93±0.03 69.67±0.86 


18 0.51 


0.63±0.10 


81.42±1.78 


14.34±0.31 


£1.00 -0.31±0.09 


16.05±0.20 


54.89±0.69 


0.15±0.03 0.30±0.19 2.53±0.03 30.77±0.40 


19 0.53 


0.48±0.09 


80.31±2.73 


14.33±0.49 


^1.00 0.01±0.19 


7.85±0.49 


23.25±1.44 


0.59±0.08 -0.27±0.09 11.84±0.18 62.42±0.97 


20 0.51 


0.32±0.69 45.91±12.77 


4.88±1.36 


1.79±0.74 ^0.70 


14.25±3.58 


6.72±1.69 


0.78±0.06 -0.24±0.09 39.85±0.16 88.41±0.36 


21 0.61 


0.30±0.22 


72.08±0.72 


12.60±0.13 


^1.00 -0.24±0.22 


19.96±0.87 


58.47±2.56 


0.53±0.25 0.39±0.21 7.96±0.69 28.93±2.50 


22 0.59 


-0.21±0.33 


78.24±4.36 


29.17±1.63 


1.28±0.54 0.53±0.16 


14.88±1.13 


27.18±2.06 


0.62±0.18-0.49±0.24 6.87±0.07 43.65±0.47 


23 0.55 


-0.40±0.21 


75.78±3.47 


26.55±1.22 


1.16±0.33 0.37±0.18 


13.19±0.90 


25.02±1.71 


0.60±0.10 -0.03±0.17 11.04±0.12 48.43±0.54 


24 0.57 


0.37±0.27 


76.08±8.91 


11.80±1.38 


1.29±0.58 0.70±0.37 


4.18±2.11 


4.20±2.12 


0.61±0.05 -0.28±0.07 19.74±0.18 84.00±0.76 


25 0.61 


0.29±0.44 


54.79±7.03 


7.64±0.98 


1.21±0.32 0.52±0.09 


19.25±1.63 


22.39±1.89 


0.79±0.07 -0.16±0.10 25.96±0.35 69.97±0.93 


26 0.68 


0.36±0.14 


90.56±1.70 


27.22±0.51 


1.23±0.21-0.56±0.30 


7.11±0.28 


38.36±1.50 


0.26±0.10-0.07±0.39 2.33±0.07 34.42±1.01 


27 0.58 


0.18±0.12 


93.24±3.59 


38.64±1.49 


1.37±0.47 ^0.70 


2.04±0.96 


4.24±1.99 


0.49±0.09-0.53±0.15 4.72±0.04 57.12±0.52 


28 0.55 


0.28±0.26 


74.91±3.04 


13.78±0.56 


1.38±0.16-0.53±0.15 


17.09±0.76 


48.99±2.18 


0.48±0.16 0.18±0.27 8.00±0.35 37.24±1.64 


29 0.54 


0.31±0.28 


84.87±8.75 


16.04±1.65 


1.86±3.84 0.67±1.74 


0.26±2.55 


0.20±2.03 


0.51±0.05 -0.25±0.10 14.87±0.07 83.75±0.40 


30 0.63 


-0.12±0.31 


60.50±1.88 


10.53±0.33 


£1.00 -0.25±0.06 


36.83±0.39 


71.34±0.76 


0.16±0.03 ^0.30 2.66±0.07 18.13±0.45 


31 0.60 


0.21±0.14 


91.49±1.07 


35.99±0.42 


1.00±0.60-0.63±0.35 


7.46±0.44 


52.18±3.05 


0.44±0.72-0.04±1.90 1.06±0.24 11.82±2.65 


32 0.86 


0.31±0.15 


81.71±0.59 


21.11±0.15 


3.91±1.84 £-2.30 


6.29±1.21 


14.02±2.69 


0.84±0.10 -0.40±0.16 11.99±0.48 64.87±2.59 


33 0.63 


0.24±0.20 


79.36±2.21 


16.37±0.46 


1.05±0.27-0.44±0.28 


12.75±0.97 


43.37±3.32 


0.55±0.26-0.09±0.52 7.89±0.57 40.26±2.88 


34 0.58 


0.27±0.15 


88.76±1.07 


26.96±0.33 


0.90±0.12-0.29±0.13 


9.97±0.22 


54.75±1.19 


0.21±0.09 0.33±0.42 1.28±0.06 18.29±0.89 


35 0.60 


0.19±0.50 


72.40±8.02 


15.09±1.67 


2.05±0.97 0.23±0.73 


13.61±2.40 


14.59±2.57 


0.53±0.08 -0.09±0.16 13.99±0.18 70.32±0.92 


c 0.22 


-0.20±0.41 


57.95±6.64 


13.05±1.50 


1.35±0.73 0.39±0.23 


22.44±1.74 


26.80±2.08 


0.72±0.09 -0.22±0.15 19.61±0.20 60.14±0.62 



measurement uncertainties. The algebraic degeneracy arises 
from number of unknowns larger than number of observ- 
ables. Unlike methods in which only a few spectral indices 
are used for synthesis, ULySS synthesize the full spectrum, 
so the algebraic degeneracy is not a problem, as the number 
of the points used in the observational spectrum far exceed 
the number of unkown parameters. Similarly, the degenera- 
cies associated with different stellar populations should b e 
relieved by fitting the full spectrum i) Jimenez et al.ll2004h . 
Furthermore, ULySS determine all the free parameters in a 
single fit so as to effectively deal with the existing degenera- 
cies between them instead of estimating the parameters in 
different steps. More expUcit explanation on th is idea of han- 
dling d egeneracies can be read in Section 2 in iKoleva et al.l 
l|2009al ). 



To check the reliability of ULySS, iKoleva et all (|2008l ) 
have analysed the spectra of Galactic clusters whose popu- 
lations are already known from colour-magnitude diagrams 
(CMD)determinations, and the two results are well consis- 
tent. So far, ULySS has been presented in several places 
to study stellar populations and its reliability assessed on 
the basis of various si mulations (jKolcva et al.,i2008 . i2009bl : 
iMichielsen et al.ll2007l '). 

To carefully verify the reliability and robustness of 
ULySS by ourselves, we analyse the "stellar populations" 
of a 'fake galaxy' made from known SSPs. Such a test has 
also been adopted by STA RLIGHT group to check th e ro- 
bustness of STARLIGHT ijcid Fernandes et all [JOOsh . In 
our test, we firstly extract three spectra with definite ages 
and metallicities from the PE model. Spectrum 1 represents 
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Table 6. Parameters from 3-population (OSP + ISP + YSP) fit with the VM model for each of the sample. This table presents the 
population histories of the 35 E+As, reconstructed from 3 episodes of star formation. Column 1 is the Galaxy Number (GN) for the 
E+A sample listed in Table 1. Column 2 is x^ value from the best fit. The following are age (in Gyr), [Fe/H] (in dcx), mass fraction 
(MF, in %), and light fraction (LF, in %) for each of the component (OSP, ISP, and YSP). 



GN x^ 


OSP 


(fixed to 12Gyr) 


ISP 


YSP 




[Fc/H] 


MF 


LF 


Age [Fc/H] MF LF 


Age [Fc/H] MF LF 




(dcx) 


(%) 


(%) 


(Gyr) (dcx) (%) (%) 


(Gyr) (dcx) (%) (%) 


1 0.40 


^0.20 


73.62±2.91 


18.68±0.74 


1.27±0.21-0.28±0.23 7.07±1.10 21.57±3.35 


0.87±0.09 0.09±0.05 19.31±0.85 59.75±2.63 


2 0.56 


0.13±0.08 


93.23±1.42 


44.21±0.67 


1.42±0.63-1.19±0.32 2.55±0.28 19.15±2.11 


0.60±0.09-0.06±0.14 4.22±0.17 36.64±1.46 


3 0.52 


^0.20 


95.19±0.71 


47.10±0.35 


0.62±0.06-0.17±0.10 4.48±0.13 43.09±1.25 


0.12±0.08 0.20±0.58 0.33±0.03 9.81±0.92 


4 0.65 


0.05±0.14 


82.89±2.61 


27.04±0.85 


^1.00 -0.05±0.15 12.97±0.67 44.02±2.28 


0.42±0.07 0.19±0.12 4.14±0.21 28.95±1.44 


5 0.45 


0.20±-0.00 


8.27±8.68 


0.89±0.93 


1.61±0.13 0.19±0.04 72.84±1.93 53.76±1.43 


0.49±0.03-0.02±0.05 18.89±0.21 45.36±0.51 


6 0.57 


0.18±0.11 


80.58±4.98 


22.88±1.41 


1.34±0.15 -0.20±0.00 10.35±1.09 30.55±3.21 


0.55±0.03 0.16±0.05 9.07±0.35 46.56±1.81 


7 0.62 


-1.13±0.18 


4.56±4.74 


1.19±1.24 


1.60±0.16 ^0.20 69.56±1.49 51.86±1.11 


0.67±0.04 -0.03±0.07 25.88±0.13 46.95±0.23 


8 0.50 


0.10±0.18 


59.46±2.58 


11.71±0.51 


1.60±0.92-0.53±0.49 6.08±1.15 11.14±2.11 


0.84±0.05 0.14±0.04 34.46±0.73 77.15±1.64 


9 0.67 


-0.95±0.24 


51.38±5.17 


19.87±2.00 


1.58±0.21 ^0.20 37.57±1.74 45.04±2.08 


0.63±0.11 -0.14±0.18 11.04±0.09 35.09±0.30 


10 0.49 


-1.09±0.20 


35.05±2.06 


18.36±1.08 


^5.00 ^0.20 55.80±2.51 16.93±0.76 


0.29±0.01 0.14±0.04 9.15±0.05 64.71±0.33 


11 0.61 


^0.20 


70.68±6.08 


15.85±1.36 


2.75±0.76 -0.37±0.14 22.21±1.66 26.56±1.99 


0.25±0.02 0.20±0.10 7.11±0.08 57.59±0.64 


12 0.52 


-0.17±0.32 33.11±10.73 


5.73±1.85 


1.49±0.09 0.14±0.07 45.80±2.90 45.84±2.90 


0.68±0.06-0.14±0.08 21.08±0.47 48.43±1.07 


13 0.70 


^0.20 


41.64±9.85 


3.50±0.83 


^5.00 -1.54±0.35 43.80±4.60 13.98±1.47 


0.12±0.00 -0.10±0.05 14.56±0.12 82.52±0.66 


14 0.64 


^0.20 


54.35±4.23 


11.68±0.91 


1.57±0.09 ^0.20 41.58±0.87 63.72±1.34 


0.38±0.03-0.04±0.12 4.07±0.07 24.60±0.44 


15 0.71 


-1.01±0.19 


45.57±3.01 


17.18±1.13 


2.17±0.26 ^0.20 29.69±1.46 21.56±1.06 


0.71±0.02 0.02±0.04 24.74±0.08 61.26±0.19 


16 0.53 


0.14±0.09 


76.61±3.46 


19.84±0.90 


1.43±0.13 ^0.20 10.47±0.87 20.11±1.66 


0.60±0.04 -0.02±0.06 12.93±0.17 60.05±0.78 


17 0.59 


0.19±0.10 


88.51±3.46 


15.11±0.59 


1.23±0.37-1.54±0.17 7.03±0.41 25.45±1.48 


0.12±0.00-0.45±0.07 4.46±0.07 59.44±0.91 


18 0.59 


^0.20 


62.20±2.47 


8.86±0.35 


^1.00 0.02±0.03 32.33±0.49 51.33±0.78 


0.16±0.02 0.14±0.12 5.46±0.06 39.80±0.44 


19 0.55 


-1.24±1.02 


7.56±4.58 


1.85±1.12 


1.64±0.09 0.19±0.06 66.84±1.53 42.87±0.98 


0.49±0.02 0.01±0.05 25.60±0.10 55.28±0.21 


20 0.51 


^0.20 


31.64±16.14 


3.47±1.77 


1.58±0.37 0.17±0.12 29.93±3.90 23.23±3.03 


0.68±0.05 -0.09±0.06 38.44±0.67 73.30±1.27 


21 0.61 


0.19±0.18 


71.42±2.10 


15.63±0.46 


0.84±0.06 0.10±0.04 27.41±0.44 75.05±1.21 


0.27±0.15-0.01±0.54 1.18±0.10 9.31±0.77 


22 0.60 


^0.20 


86.07±3.29 


29.72±1.14 


1.24±0.15 -0.20±0.10 11.63±0.57 46.23±2.26 


0.32±0.07 0.19±0.29 2.30±0.11 24.05±1.14 


23 0.53 


-0.52±0.24 


52.46±6.69 


14.75±1.88 


1.44±0.09 0.10±-0.00 35.19±1.99 46.17±2.62 


0.51±0.05 0.17±0.05 12.35±0.24 39.08±0.76 


24 0.60 


^0.20 


75.12±2.54 


15.15±0.51 


2.53±0.95-1.29±0.37 7.78±0.68 15.30±1.33 


0.54±0.03 0.01±0.05 17.10±0.20 69.55±0.82 


25 0.61 


0.17±0.14 


70.58±2.69 


15.57±0.59 


0.92±0.07 0.17±0.08 20.91±0.87 49.71±2.06 


0.59±0.08 O.OlitO.lO 8.51±0.36 34.72±1.49 


26 0.77 


0.16±0.07 


91.24±1.32 


32.26±0.47 


1.23±0.13-1.51±0.28 3.88±0.19 28.16±1.39 


0.42±0.02 0.20±-0.00 4.88±0.12 39.58±0.94 


27 0.64 


^0.20 


92.86±0.91 


37.35±0.37 


^1.00 -0.03±0.18 2.42±0.26 11.23±1.23 


0.55±0.06-0.72±0.08 4.72±0.08 51.42±0.89 


28 0.59 


-0.66±0.19 


53.84±4.23 


16.91±1.33 


1.53±0.11 ^0.20 31.22±1.36 36.58±1.59 


0.56±0.04 0.01±0.09 14.95±0.11 46.51±0.33 


29 0.54 


-0.54±0.23 


58.11±9.55 


13.66±2.24 


1.61±0.35 ^0.20 25.81±2.87 23.15±2.58 


0.32±0.01 ^0.20 16.07±0.09 63.19±0.37 


30 0.62 


-1.03±0.13 


58.69±2.43 


19.72±0.82 


1.24±0.05 0.03±0.06 32.93±0.82 44.27±1.10 


0.31±0.01 ^0.20 8.39±0.07 36.01±0.31 


31 0.66 


-0.47±0.14 


70.69±4.39 


28.42±1.76 


1.59±0.21 0.14±0.08 24.87±1.34 41.16±2.22 


0.32±0.02 ^0.20 4.44±0.07 30.42±0.48 


32 0.86 


-0.02±0.13 


85.35±1.81 


30.80±0.65 


2.61±0.66 ^-1.68 6.44±0.57 20.54±1.81 


0.57±0.07 O.OOitO.13 8.20±0.20 48.65±1.18 


33 0.65 


0.19±0.17 


77.35±4.17 


18.26±0.99 


1.03±0.21 -0.22±0.13 14.95±1.24 44.67±3.70 


0.52±0.08 0.08±0.10 7.70±0.57 37.08±2.73 


34 0.57 


^0.20 


86.94±1.63 


29.53±0.55 


0.93±0.08 -0.05±0.13 10.57±0.35 45.15±1.51 


0.32±0.04 0.20±0.10 2.48±0.10 25.32±0.97 


35 0.61 


0.09±0.14 


78.63±3.07 


19.70±0.77 


0.97±0.10 0.02±0.11 15.92±0.70 42.08±1.86 


0.32±0.02 0.20±0.10 5.45±0.16 38.22±1.11 



an old stellar population with age=10 Gyr and [Fe/H]=- 
1.0 dex. Spectrum 2 represents an intermediate-age stellar 
population with age=2500 Myr and [Fe/H]=-0.5 dex. Spec- 
trum 3 represents a young stellar population with age=500 
Myr and [Fe/H]=0.2 dex. For a convenient use, we represent 
Spetrum 1, 2 and 3 as Old, Inter and Young, respectively. 
Secondly, we generate a composite 'fake galaxy' spectrum by 
the linear combination of O.lxyoung, 0.3xinter and 0.6xold 
(composite spectrum=0.1xyoung+0.3xinter-|-0.6xold). Fi- 
nally, we use ULySS to synthesize the 'fake galaxy' and de- 
rive its components. The results have a good agreement with 
the known components (See Table 8). 

In addition, we make 500 Monte-Carlo sinmlations 
which could visualize the degeneracies for this solution. 
In every simulation, we add a random Gaussian noise to 
the 'fake galaxy' and then synthesize it by ULySS. From 
Monte-Carlo simulations presented in Figure 6, the YSPs 



Table 7. Mean light fraction, mass fraction and age for the three 
populations of the 35 brightest E-l-A galaxies, illustrated respec- 
tively from the PE and VM model. 



Model 



PE 

VM 



Mean light fraction 
YSP ISP OSP 



51% 30% 
46% 35% 



19% 
19% 



Mean mass fraction 
YSP ISP OSP 



11% 14% 
12% 25% 



75% 
63% 



Mean age (Gyr) 
YSP ISP OSP 



0.54 1.82 12 
0.46 1.82 12 



(the left blue gathering area) concentrate at its true place 
(age=500Myr and [Fe/H]=0.2dex), compactly. However, the 
ISPs (the middle green gathering area) and OSPs (the right 
red gathering area) scatter a little, but anyway they can be 
divided into two different components. The gathering cen- 
ters (plus sign) of simulations are consistent with our fit 
results, which are the closest to the true value (see Table 8 
for comparison). Thus we believe that ULySS can at least 
give the most likely solutions. This is enough for us, because 
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Table 8. Comparison of the known and synthesized information 
for 'fake galaxy'. 





The known 'fake 


galaxy' 




SSP 


age(Gyr) 


[Fe/H](dex) 


Weight 




Old 


10.00 


-1.00 


0.60 




Inter 


2.50 


-0.50 


0.30 




Young 


0.50 


0.20 


0.10 






The 


synthesis results by ULySS 




SSP 


age(Gyr) 


[Fe/H](dcx) 


Weight 


light fraction 


Old 


10.00±1.77 


-1.00±0.12 


0.60±0.03 


28.64±1.30 


Inter 


2.50±0.32 


-0.50±0.11 


0.30±0.01 


34.33±1.40 


Young 


0.50±0.02 


0.20±0.03 


0.09±0.00 
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Figure 6. 500 Monte-Carlo simulations for the 'fake galaxy'. 
From top to bottom, in each picture, Young components are 
represented by blue small squares (the left gathering area). 
Intermediate-age components are represented by green small 
squares (the middle gathering area). Old components are rep- 
resented by red small squares (the right gathering area). 



we just would like to find out the dominant population and 
some quanlitative conculsions. 



4.2 Validity of our solutions 

To check the reliability of our solutions, we deal the 35 E-l-A 
sample galaxies as a whole, and combine their spectra to be 
one composite spectrum by using IRAF task-SCOMBINE. 
This combined spectrum could be a good representative of 
E-l-A galaxies. 



For this combined spectrum, we do the same 3- 
population fit with the PE model as made in §3.3. The re- 
sults are listed in Table 5 (the last line). Then we make 
500 Monte-Carlo simulations for this combined spectrum to 
verify our solutions through visualizing the degeneracies (see 
Figure 7). During each simulation, a different random Gaus- 
sian noise is added to the galaxy and then the 3-population 
fit with the PE model is made. In Figure 7, the 500 results 
of each component scatter, but they all respectively gather 
towards a center (black plus sign) which is consistent with 
our fit results. The three components (young, intermediate 
and old) represented by different colours (blue, green and 
red) and gathered respectively in the left, middle, and right 
area in every picture can be clearly seperated . Although 
scatters in each component exists and then the age of each 
component can move around within the gathering area filled 
with its colour, in the top picture (relation of age and light 
fraction), the distribution of YSPs (the left blue region) are 
statistically at the top, ISPs (the middle green region) at 
the middle, and OSPs (the right red region) at the bot- 
tom, along the vertical axis (light fraction). This means in 
a statistical sense, small changes in age do not affect the 
dominance relationship in light among YSP, ISP and OSP. 
Similarly, the middle picture convinces that the OSPs sta- 
tistically contribute the most to the mass, and the YSPs 
contribute the least. These solutions from simulations are 
consistent with our previous solutions for the 35 brightest 
E-l-A galaxies, that is, to light, YSPs are dominance, then 
ISPs and finally the OSPs; to mass, OSPs are dominance, 
then ISPs and finally the OSPs. 

For a more clear comparison, we list the fit result (F) 
copied from the last line of Table 5 and the mean result from 
our 500 simulations (S) for this combined spectrum in Table 
9. From Table 9, we can see that for the combined spectrum, 
the mean results from the simulations are generally consis- 
tent with the fit results listed in Table 5 (the last line). For 
this combined spectrum, no matter which result is consid- 
ered, F or S, YSPs contribute about 60%, ISPs about 30% 
and OSPs about the rest 10% of the light. However, OSPs 
contribute about 50%, ISPs about 30% and YSPs about the 
rest 20% of the mass. This is completely corresponding with 
our previous statistical solutions. Such consistency can un- 
derpin the correctness of our solutions, since this analysed 
combined spectrum represent all the E+A sample galaxies 
as a whole. 

Another way to check the reliability of our solutions is 
to fix the fitting technique but change the population models 
to fit our sample. In our work, the PE and VM model are two 
different and independent population models. We have fitted 
all of our sample using both of the models. The agreement 
of the results from the two models is satisfactory (see details 
in § 3.2 and § 3.3). 

Comparing our solutions with those from other work 
is also a good way to verify the reliability of our solutions. 
Our work has shown in Table 7 that the optical light of 
E+As are dominated by YSPs (<lGyr), the average light 
fraction of which is the most of all (51% illustrated with the 
PE model and 46% illustrated with the VM model). The 
average light fraction is less in the ISPs (30% with the PE 
model and 35% with the VM model), and the least in OSPs 
(19% with the PE model and 19% wit h the VM model) 
Our solutions are consistent with those of lPracv et al.l (|2009l ) 



12 W. Du et al. 

Table 9. Comparison between the fit results (F) and the 500 simulations results (S) for the combined spectrum. The first line shows 
the optimal results copied from Table 5, and the second line shows the mean results (the position of the plus sign in Figure 7) from 500 
simulations. 
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Figure 7. 500 Monte-Carlo simulations for the combined spec- 
trum. From top to bottom, in each picture, Young components 
are represented by blue small squares (the left gathering area). 
Intermediate-age components are represented by green small 
squares (the middle gathering area). Old components are rep- 
resented by red small squares (the right gathering area). 



for the stellar population study of the 2dF E+A sample 
galaxies. Their work has revealed that the composite model 
spectra that best fit the E-l-A spectra are dominated by YSP 
templates (<lGyr). The agreement between the solutions 
from our work and other work on E+As further underpins 
the reliability of our solutions. 

All this series of tests and comparisons could verify that 
our solutions are robust. 



5 RECONSTRUCTION OF THE STAR 
FORMATION HISTORIES OF THE 
SAMPLE SFR 

The fundamental information that characterizes CSPs is its 
SFH, that is the evolution with time of the amount (i.e. total 
mass) of stars formed (SFR). The traditional and simplest 
way to probe the galaxy evolution is to measure the observed 
SFR in galaxies. 

ULySS does not imediately give the SFR after fitting, 
but it gives the optimal weights of the different SSP com- 
ponents. To derive the SFR we have to assume two points. 
Firstly, the individual SSPs approximate the SFH over a pe- 
riod of time. Secondly, we regard the SFR as a constant in 
each time box. If only these two assumptions are taken into 
account, we can then calculate the SFR. 

First of all, we have to choose the limits of these time 
boxes. We set limits at medium location between the loga- 
rithm ages of the components, and derive the outer limits 
so that the age of the extreme components sit in the mid- 
dle of their boxes. This approach is the smoothing strategy, 
which associates a duration to each burst and then smooths 
the instantaneous burst in its assuming duration time box 
by convolving with a Gaussian function. The choices of the 
time box above is nearly the centers of the Gaussian func- 
tion. We consider the 12Gyr in the past as the beginning 
of the star formation. Such choices for tim e bo xes have al- 
ready been adopted in lKoleva et al.l (|2009br ) and lAsari et al.l 
(! 20071) to get generally smoothing star formation histories. 

Then, the smoothing SFR in each time box is calculated 
out as the weight of the correlated stellar component divided 
by the length of the time box. Thus, we obtain the smooth- 
ing SFRs for all of the 35 E-f As. In Figure 8, we present 
the SFRs for each of the sample as a function of evolution- 
ary time. For each of the sample, the SFR illustrated from 
the fit with the PE model is represented by the solid line 
while SFR illustrated from the VM model is represented by 
the dashed line. When we see Figure 8, all the assumptions 
above have to be taken into account. This figure at least al- 
lows a direct comparison of the histories reconstructed with 
the two different models. 

In Figure 8, the SFR for each of the sample have good 
agreement between the PE model and the VM model. Al- 
most all of the sample galaxies have the relatively low SFR 
in the intermediate evolutionary stage, and less than 0.1 so- 
lar mass of stars is formed stably every year in this stage. 
For 11 E-f A galaxies (Galaxy number 2, 3, 10, 11, 13, 17, 
26, 27, 31, 32 and 34), the early stage is important for their 
formation, as they have high SFRs, which are greater than 
0.1 solar mass every year in the early period. For other 11 
E+A galaxies (Galaxy number 1, 4, 5, 6, 8, 9, 12, 15, 16, 
19 and 20), the SFRs during the late evolutionary stage 
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Figure 8. SFR illustrated from the fit with the PE model (in solid line) and the VM model (in dashed line) for eaeh of the sample. 
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are quite high, and at least 0.15 solar mass of stars can be 
formed in this stage. Such a phenomenon could be due to 
the starburst in the E+A galaxy which takes place recently 
and is truncated in a short time. For these 11 E+As, it is 
probable that the timescale for the lateset formation phase 
is much shorter than assumed and the instantaneous SFR 
of this phase is much higher. These are just our explanation 
and we expect to explore more in our future work. 



6 EVOLUTIONARY SCENARIO 

It is commonly accepted that E+A galaxies are the trans- 
formation phase from late-type galaxies to early-type galax- 
ies, and they are likely to evolve in to E/SO galaxies, which 
are the typical ea rly-type galaxies (jYang et al.l |2004| . 120081 : 
iFalkenberg et al.ll2009i V However, the true nature and the 
validity of the link between E-fAs and E/SOs remain un- 
certain, while some auth ors make their wo rks to increase 
the reliability of this link. lYang et al.l ((2003) have used high 
spatial resolution images to probe the detailed morpholo- 
gies of the local E-l-A galaxies, and they find that the prop- 
erties of morphologies, color profiles, scaling relations, and 
young star clusters of E+As are either consistent with those 
of E/SOs, or indicate that E-fA galaxies are caught in the 
act of transformin g from late-type to early-type galaxies. 
IPracv et al.l ((20091) also prove the evolutionary link from the 
evidence of the morphology. Their E+A galaxies are con- 
sistent with being early-type systems based both on radial 
surface brightness profiles and visual morphological classifi- 
cation. Besides, the ubiquitous rotation found in their E+A 
galaxies provides that E+A galaxies are completely consis- 
tent with the 'fast rotator' population of early-type galaxies. 
This gives another evidence to prove the evolutionary link 
between the E+As and E/SO (early- type galaxies). 

In this work, we expect to give another evidence from 
ages of the stellar populations to support that evolution- 
ary link between E+As and early-type galaxies. If there is 
sure to exist the evolutionary link from E+As to early-types 
(E/SOs), we would expect the early- types have, in general, 
older old stellar population ages than E+As. 

To complete it, we plan to fit our E+As and 34 
SPSS early-type g a laxies (ellipticals) randomly selected 



from [Rogers et al.l ([20091) with 2-population models by 
ULySS since we hope to seperate their components into two 
groups in age (one is the old age group, and the other is 
the young age group), and then we compare the equivalent 
age of the old components from our E+A sample and the 
selected ellipticals. It is expected that the equivalent age of 
old components in ellipticals is older than that in E+As. If 
this hypothesis is true, an evidence in stellar population ages 
can be obtained to support the evolutionary link. While we 
are fitting the E+A sample, we choose SSPs aged between 
lOMyr and 2Gyr as young components, and SSPs aged be- 
tween 2Gyr and 16Gyr as old components, and we fit each of 
the E+A sample with these two components by using ULySS 
with the PE model. The age limits for two components are 
fine for all of E+A galaxies. We obtain the equivalent ages 
of old component and young component for each of the sam- 
ple, and the validity of the solution is under the robustness 
of ULySS (§ 4.1) and is also tested by the same Monte-Carlo 
simulations as described in § 4.2. For elliptical galaxies, if 
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Figure 9. Distribution histogram of old stellar population ages 
for E+As and ellipticals from 2-population fit with the PE model. 
We use bin=lGyr for this histogram. The solid line shows the age 
distribution of old components for E+As, and the dashed line 
shows for ellipticals. The mean old component age for E+As is 
T.OGyr, and for ellipticals is ll.SGyr. 



we use the same age limits for components, the age of young 
component would always reach the upper limit (2Gyr) after 
fitting, which reveals an older than 2Gyr age is required for 
young component of ellipticals. So we choose SSPs aged be- 
tween lOMyr and 5Gyr as young components and SSPs aged 
between 5Gyr and 16Gyr as old components to fit ellipticals. 
These new limits are fine for all the ellipticals. 

We obtain the equivalent age of the old component (to) 
for each of the E+A galaxies and elliptical galaxies. Our 
results show that to of 30 E+As (30/35=86%) distribute 
lower than lOGyr, and the average age of the old stellar 
populationsof E+As is concentrating on 7.0Gyr. However, 
for ellipticals 30 of them (30/34=88%) have to larger than 
9Gyr, and the average age of the old stellar populations of 
the ellipticals is 11.5Gyr, which is much older than that of 
E+As (see distributions in Figure 9). 

Therefore, our workis consistent with the proposed evo- 
lutionary link from E+As to early-types (E/SOs) in the as- 
pect of stellar population age. 



7 SUMMARY AND CONCLUSIONS 

In this paper we have presented a detailed stellar popula- 
tion synthesis work on our 35 brightest E+A galaxies se- 
lected from Goto's E+A galaxies catalogue from SDSS DR5 
(|Adelman-McCarthv et all \200i\ ). These E+As have pet- 
rosian magnitudes in r-band less than 16 mag and spectro- 
scopic S/N ratio in the continuum of the r-band wavelength 
range greater than 30. The redshifts of them are all within 
z = 0.14. Considering our selection criteria (brightest), the 
selected 35 E+A sample are good representatives of the full 
564 E+As from Goto's catalogue in terms of the distribu- 
tions of redshift, absolute magnitude, diameter, velocity dis- 
persion, SSP-equivalent age and metallicity. 

We have fitted the 35 E+A sample with a SSP model 
and CSPs model respectively from the PE and the VM 
model. The SSP-equivalent age, metallicity and velocity dis- 
persion of each of the sample can be obtained by a SSP fit. 
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which has a mean SSP-equivalent age of 0.87Gyr with the 
PE model and l.UGyr with the VM model, a mean [Fe/H] 
of -0.04 dex and -0.38 dex, and a mean velocity dispersion 
of 159 km/s and 141 km/s (see Figure 2). The results from 
the two different models are generally consistant. 

For the CSP fit, we use the combination of YSP + ISP 
-|- OSP and get the optimal parameters of YSP, ISP and 
OSP for each of the sample. We have detected the 12Gyr 
OSPs with insignificant light fraction for all of the sample 
galaxies. However, the YSP components (< IGyr) dominate 
the E-|-As in light, but contribute very little to the mass of 
E-l-A galaxies. Our conclusions are that for the E+A sample, 
YSPs are statistically the dominance to the light, then ISPs 
and finally OSPs. However, to mass, OSPs are statistically 
the dominance, then ISPs and finally YSPs. 

We verified the reliability and robustness of our solu- 
tions from 4 aspects: the robustness of the technique we 
used, the validity of our solutions, the consistency of the 
results from two independent population models and the 
agreement with other studies. 

We have reconstructed the SFH for each of the E-l-A 
sample by giving their smoothing SFR. For 11 E-|-As, the 
early stage is the significant period for star formation. How- 
ever, for another 11 E-|-As, their SFRs in the late evolu- 
tionary stage are very high, which possibly arises from the 
extremely sudden truncation of the recent starburst. 

At last, we additionally compared the stellar popula- 
tions of the E-|-As with those of the early-type gala:xies. We 
fitted our 35 E-|-As and the 34 early-types randomly selected 
from SDSS with 2-population model (OSP + YSP) from the 
PE model. The results show that the ellipticals have in gen- 
eral older OSP ages than E-f As. This gives at least another 
evidence in the aspect of stellar population age to support 
the proposed evolutionary senario of E-|-As evolving into 
early- types (E/SOs). 
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